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Abstract 

Pioneering atmospheric muon neutrino experiments have demonstrated the near-maximal 
magnitude of the flavor mixing angle 622,- But the precise value of the deviation D = 1/2 — 
sin 2 #23 from maximality (if nonzero) needs to be known, being of great interest - especially 
to builders of neutrino mass and mixing models. We quantitatively investigate in a three 
generation framework the feasibility of determining D in a statistically significant manner 
from studies of the atmospheric Vn, survival probability including both vacuum oscillations 
and matter effects. We show how this determination will be sharpened by considering the 
up-down ratios of observed v^- and ^-induced events and the differences of these ratios 
in specified energy and zenith angle bins. We consider 1 Megaton year of exposure to a 
magnetized iron calorimeter such as the proposed INO detector ICAL, taking into account 
both energy and zenith angle resolution functions. The sensitivity of such an exposure and the 
dependence of the determination of D on the concerned oscillation parameters are discussed 
in detail. The vital use of matter effects in fixing the octant of 623 is highlighted. 
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1 Introduction 



Neutrino oscillation studies have come of age after the recent results from the Super-Kamiokande 
(SK) [1], K2K [2], SNO [3] and KamLAND [4] experiments. The existence of oscillations in atmo- 
spheric [5] and accelerator generated [2] u^s and P^s with a modulus of squared mass difference 1 
|Am^| ~ 2.1 x 10~ 3 eV 2 (1.4 x 10" 3 eV 2 < |Am^| < 3.3 x 1(T 3 eV 2 at 3a) and a mixing an- 
gle 023 ~ 45° (35° < 9 23 < 54° at 3a) as well as in solar [3, 6] (reactor [4]) z/ e 's (P e 's) with a 
squared mass difference 1 Am 2 21 ~ 8 x 10~ 5 eV 2 (7.4 x 10~ 5 eV 2 < Am 2 21 < 9.5 x 1(T 5 eV 2 at 3a) 
and a mixing angle 9 12 ~ 34° (30° < 9i 2 < 40° at 3a) is now an accepted fact [7]. The goals 
of the next generation of experiments include both increased precision in the already measured 
oscillation parameters and the determination of the still unknown quantities such as 13 (presently 
bounded by < #i 3 < 12° at 3a [7, 8]) as well as the sign of Am| x and eventually the CP violating 
phase S in the PMNS matrix. An issue of great importance right now is the precise value of 9 23 
which dominantly controls v jX -v T oscillations. The presently allowed range of #23 does not enable 
one to distinguish it from its maximal 2 value 7r/4. A naturally arising question is: how can one 
significantly narrow down this range in the future and possibly detect a nonzero deviation 

£>=I_ sin 2 23 (!) 

from maximal mixing. This is what we consider in the context of forthcoming 3 studies of atmo- 
spheric z/^'s and P^s. 

Let us mention some current theoretical ideas regarding 9 23 . A simple way of understanding 
a maximal value for 9 23 is to invoke the concept of a /i-t exchange symmetry [10]. Though 
this immediately leads to #23 = 7r /4, it also yields at the same time a vanishing 9 13 and CP 
conservation in the neutrino sector. While the last two features could in principle turn out to be 
the facts of Nature, our current hopes are to the contrary. A nonzero and measurable #13 as well 
as the occurrence of CP violation in the neutrino sector can lead to much that is interestingly 
new in neutrino physics and can be studied in laboratory experiments with reactor, accelerator, 
atmospheric, solar and supernova neutrinos. The general expectation [10], therefore, is that any ji- 
t interchange symmetry, if present, must be broken so as to allow a nonvanishing 13 and (possibly 
observable) CP violation in the neutrino sector. Such a breaking would generally cause D to be 
different from zero. However, its magnitude and sign would be given by the yet undetermined 
symmetry breaking mechanism. Another idea has been that of quark- lepton complementarity [1 1] . 
In this approach a bimaximal neutrino Majorana mass matrix is "rotated" by a unitary matrix 
diagonalizing charged left chiral leptons in generation space. The latter is postulated to be the 
same as the CKM matrix resulting in a reduction of 9\ 2 from 7r/4 by an amount comparable to 
the Cabibbo angle 9c ~ 13° in conformity with observation. In this case #23 also is expected to 

*We work in the conventional picture of three neutrino flavors with CPT conservation assumed. We also employ 
standard notation where mfj — mf — mj and Oij plus S are defined by the PMNS matrix of our Eq. (2) in §2. 

2 Two flavor v^jvy, oscillations in vacuum are controlled by #23 only through the factor sin 2 26*23 which is high- 
est when 023 has the maximal value 7r/4. Matter effects and Am^-driven subdominant atmospheric neutrino 
oscillations have different dependence on #23. 

3 An analysis of the allowed range of D, based on extant atmospheric + data, is available in [9]. 
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be reduced from 45° by an amount of the order of sin - S23 ~ 2.4° so that D is expected to be 
positive and in the vicinity of 0.04. However, variants [11] of this hypothesis exist with D predicted 
to be significantly different from 0.04. There are many other model predictions, utilizing flavor 
symmetries, GUT relations and textures, covering \D\ from 0.005 to > 0.16, as listed in Table II 
of [12]. 

We thus see that an experimental measurement of D would be of great interest to the builders 
of neutrino mass matrix models. Apart from the magnitude of |D|, sgn(D) is also of importance. 
Its determination would fix the octant of #23, i-e- whether #23 exceeds or is less than ir/4. With the 
goal of determining D, one can explore different experimental options. A detailed survey of the 
capabilities of forthcoming and futuristic accelerator based long baseline neutrino/antineutrino 
experiments along this line has already been given in [12, 13]. Here we want to consider future 
atmospheric neutrino studies guided by our knowledge of extant ones [9]. A major new venture 
in this direction will be the proposed Hyper-Kamiokande/UNO/MEMPHYS type of experiments 
[14] with Megaton water Cerenkov detectors. Detailed studies [15, 16, 17] have been made (among 
other things) of the + survival probability, information on which will emanate from such a 
detector. While this kind of data can yield information also on D [9], it will be largely from the 
vacuum oscillations of and and will be less sensitive to matter effects. Depending on the 
mass ordering of the neutrinos (normal or inverted), larger matter effects appear in the neutrino or 
antineutrino channel for long baseline distances (L > 1000 km). In addition, the sign of the earth 
matter effect term in the difference of survival probabilities AP W (as defined in §2) depends on the 
energy and zenith angle of the atmospheric v^jv^. Since water Cerenkov atmospheric neutrino 
experiments measure the + survival probability and since their poor detector resolution 
means that the data collected would have to be grouped into very broad sub-GeV and multi-GeV 
energy bins as in the SK data sample, matter effects in this class of experiments will be partially 
washed out. In order to maximize the extraction of information from the individual survival 
probabilities of atmospheric u^s and P M 's, propagating both through vacuum and through earth, 
one needs to measure them separately as a function of energy and baseline length. In fact, the 
individual survival probabilities will be measurable in a large magnetized iron calorimeter such as 
the ICAL detector of the proposed INO experiment [18]. This experiment is expected to have a 
good detector resolution allowing fine binning of the data in energy and zenith angle, essential for 
gleaning out the D dependence of earth matter effects. Our aim is to focus on the feasibility of 
utilizing the data in future from such a detector to determine D. 

There have already been several studies [19, 20, 21] of the possibility of determining the nor- 
mal or inverted nature of the mass ordering of neutrinos, i.e. S(7n(Am| 1 ), in magnetized iron 
calorimetric detectors. The question of the determination of D has been touched upon [20] but 
not analyzed in detail. In this paper we carry out a detailed study of the feasibility of a mea- 
surement of both the magnitude and the sign of D from data simulated for 1 Megaton year of 
exposure of a large magnetized iron calorimeter, such as ICAL, to atmospheric muon neutrinos 
and antineutrinos. In principle, a magnetized iron calorimetric detector in an underground ob- 
servatory (such as the proposed INO laboratory) could also study other types of events: upward 
going muons and electron neutrino events. As a first step, we would like to restrict our studies 
to fully contained events only and hence do not consider upward going muons, which are much 
harder to incorporate. Furthermore, the large thickness [18] of the iron plates of the proposed 
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ICAL detector precludes the detection of electrons and will not trigger on to events induced by 
u e s and P e s. This is why we confine ourselves to atmospheric v^v^ studies. Furthermore, ICAL 
is designed to have a detection threshold of neutrino energy = 1 GeV and therefore cannot study 
sub-GeV events. This is why we consider the multi-GeV regime only. We take a three generation 
system and use exact numerical solutions of the neutrino (and antineutrino) equations of motion in 
the atmosphere and in earth matter using the PREM [22] density profile as an input. We assume 
what we regard as reasonable energy and zenith angle resolution functions for the detector and 
concentrate on up-down asymmetry ratios [20] Un/D n and Ua/Da for v^- and z/^-induced events 
in appropriate energy and zenith angle bins. We are then able to pinpoint the precise bins from 
where the utilizable information on both vacuum oscillations and matter effects can be extracted - 
leading to a determination of D. We also show how the use of matter effects through the difference 
Un/Dn — Ua/Da will enable one to fix the sign of D and resolve the octant ambiguity in # 2 3- 

The paper is organized in the following way. Section 2 contains a slightly simplified and 
approximate analytical treatment of the survival probabilities of is^s and z/ M 's propagating in 
matter of constant density; this is to bring out the basic physics issues in question, specifically, 
highlighting the phenomenon of resonance. In §3 we present our numerical results on the survival 
probability (after propagation in earth matter) as it varies with the energy E and baseline length 
L; in particular, the occurrence of extrema in this variation is highlighted through the definitions 
of SPMIN1, SPMAX and SPMIN2 which appear. Section 4 comprises the methodology for the 
extraction of the up-down asymmetry ratios Un/Dn and Ua/Da from the simulated data. The 
details of our statistical (x 2 ) analysis are presented in §5. The determination of \D\ from the 
simulated data and the estimated intervals on both sides of D — for which this will be possible 
at the 3cr level are given in §6. The procedure for the fixation of the sign of D (and hence the octant 
of #23) is described in §7. In §8, which contains our summary and conclusions, we also provide 
a discussion of the sensitivity of an ICAL-like detector to D, as compared to the forthcoming 
long baseline and water Cerenkov detector studies with accelerator generated and atmospheric 
neutrinos respectively. 

2 Approximate analytical treatment of v^/ survival prob- 
ability in matter 

In the standard parametrization, the Pontecorvo-Maki-Nakagawa-Sakata [23] matrix is given by 



In Eq. (2) Cjj = cos^-, Sij = sin^- and, while the phase S has been retained, the two Majorana 
phases have been ignored since they do not contribute to neutrino oscillations. The effective 
Hamiltonian for such oscillations in matter with varying density can be expressed in the flavor 
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basis as a function of the path length x of the neutrino from its source. Thus 



H(x) = —U | Am 
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In Eq. (3), Am 2 , 



mj 



m 2 , rrii and E being the mass of the ith (physical) neutrino and the 



neutrino energy respectively, whereas the potential V(x) is given by 

V{x) = V2G F N e (x). (4) 
Here N e (x) is the electron density of matter, so that the potential can be rewritten as 



V(x) = (7.56 x 1(T 14 ) 



P(s) 
gms/cc 



Y e (x) eV, 



p(x) being the density of earth matter in the path of the neutrino and Y e (x) 
number of electrons per nucleon in the same. 

The neutrino evolution operator S(x, 0) has the matrix element 

S lP (x,0) = (i/ 7 (x)|^(0)), 

P, 7 being flavor indices. This obeys the evolution Equation 

.dS^(x, 0) 



dx 



= [H(x),S(x,0)] l/3 . 



(5) 

0.5) being the 
(6) 
(7) 



The probability for a neutrino flavor transition (3 — * 7 at a baseline length L is given by 

P Pl (L) = P[u p (0) - i/ 7 (L)] = \S^(L,0)\ 2 . (8) 

In case the earth matter density p is taken 4 to be constant between the production and detection 
points of the neutrino, the potential V also becomes a constant. It is then useful to diagonalize 
the Hamiltonian of Eq. (3) with eigenvalues Ai 5 2,3(2-E) _1 . 



H=—U M 
2E 



/Ai 
A 2 I V 

V A 3 



Aft 



(9) 



U being the lepton mixing matrix in matter. Then one can write 



(10) 



i=l 



A m 2 r 



2£ 



5 



E in GeV 


L in km 


A 


AA 


1 


1000 


0.099 


0.252 


1 


3000 


0.125 


0.953 


1 


5000 


0.131 


1.673 


1 


7000 


0.155 


2.754 


1 


9000 


0.171 


3.908 


1 


11000 


0.229 


6.409 


3 


1000 


0.297 


0.252 


3 


3000 


0.375 


0.953 


3 


5000 


0.395 


1.673 


3 


7000 


0.465 


2.754 


3 


9000 


0.513 


3.908 


3 


11000 


0.688 


6.409 


5 


1000 


0.495 


0.252 


5 


3000 


0.626 


0.953 


5 


5000 


0.659 


1.673 


5 


7000 


0.775 


2.754 


5 


9000 


0.855 


3.908 


5 


11000 


1.147 


6.409 


7 


1000 


0.693 


0.252 


7 


3000 


0.876 


0.953 


7 


5000 


0.922 


1.673 


7 


7000 


1.084 


2.754 


7 


9000 


1.197 


3.908 


7 


11000 


1.606 


6.409 


9 


1000 


0.892 


0.252 


9 


3000 


1.126 


0.953 


9 


5000 


1.185 


1.673 


9 


7000 


1.394 


2.754 


9 


9000 


1.538 


3.908 


9 


11000 


2.064 


6.409 



Table 1: Variation of A and A A with E and L assuming p to be constant and equal to 4.52 
gms/cc. 
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More specifically, the muon neutrino and antineutrino survival probabilities are given by 



^,2,^,2^2^1 \2 L + |f/ M |2|f/ M|2 sin 2 



T 1 Mil i w m3I ° 



+ l^| 2 |^| 2 sm 2 ^^Ll, (12) 



P m (L) = P^V ^ -V). (13) 

The elements of U M , appearing in Eqs. (9), can be parametrized in the same way as the 'vacuum' 
mixing matrix U of Eq. (2) except that one needs to use the corresponding values of the angles 
and the phase in matter, i.e. 0^ — > Of- and 5 — > 5 M . 

Tractable analytic expressions P a p(L) emerge only after some additional approximations. For 
the purpose of displaying the relevant analytic expressions, in addition to assuming the constancy 4 
of p and hence of V, we choose to neglect 0(sf 3 ) and 0(a 2 ) terms, a being Am^/Am^. From 
what was discussed in §1, we already know that 

4 < 0.008, (14) 
a 2 ~ 0.001. (15) 

Therefore these approximations do not generally make any significant practical difference from 
the exact numerical results, as will be seen later. On the other hand, they do enable us to display 
the two quantities of our interest, namely P^ c , which is the muon neutrino survival probability 
in vacuum, and AP MM (L) = P MM (L) — P^(L), which is the difference between the muon neutrino 
and antineutrino survival probabilities in matter, in analytic form, making their physical features 
rather transparent. The eigenvalues of A ii2)3 and the elements of the matrix U M of Eq. (9) can 
be calculated within the above mentioned approximations following the method described in Ref. 
[24]. They can then be substituted in Eq. (11) to compute P^ C (L) and AP^(L) for atmospheric 
neutrinos. 

In order to display the desired analytic expressions, we first define a dimensionless quantity 

^ 4E~' (16) 
The expression for P^ c can now be given as [24] 

^T^ 1 -444(1- 44) sin 2 A 

+ 4aci 2 C2 3 (ci2C23 - 2si 3 si 2 4 cos 5) A sin 2A + 0(a 2 , sf 3 ). (17) 



4 This is an inaccurate assumption for the passage of GeV and multi-GeV atmospheric neutrinos through the 
earth when L exceeds 1000 km. However, we need this assumption only for displaying analytic expressions. Our 
numerical results are obtained without assuming a constant p, and actually, with the PREM [22] earth matter 
density profile as an input. Moreover, they are correct to all orders in a = Am^/Ara^ and S13. 
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For the description of matter effects, it is convenient to define another dimensionless quantity 



A = 



2EV 
Am 2 



(18) 



31 



For atmospheric muon neutrinos of GeV and multi-GeV energies and low to large pathlengths, the 
magnitude of A varies from about 0.1 to about 2.1. On the other hand, the energy independent 
AA varies from about 0.25 to about 6.4 for p ~ 4.52 gm/cc, cf. Table 1, where the values of these 
parameters have been given for various choices of E and L. The expression for AP^(L) now reads 
[24] 



4s 2 s 2 

^*13*23 



sin 2 d 



A) A sin 2 (l - A) A 



(i + Ay 



[l-Af 



q 2 2 2 
-OS 13 C 23 S 2 3 



. , . fsin(l + A)A sin(l-A)A 
sin A cos AA l —± -i ^ — ^~ 



(1 + ^4) 



i - Ay 



+ A sin 2 A 



A 



I- A 2 



-8aS 13 C 12 S 12 C23S23 cos 5 



cos A 



sin AA \ sin(l + A) A sin(l - A) A 



A 



1 + A 



I- A 



2 A sin AA 
i 



2 \ A I A Oll-L 1 ■ 

■ S 2 3 ) Sm A { S1IlA 1 ~ A 2 + 



+0(a 2 ,4). 



'cos(l + A)A cos(l-A)A N 
v 1+A I- A , 



(19) 



This expression vanishes for S13 = 0, clearly showing the need for the latter to have a nonzero value 
in order to make AP^ nonzero. If the limit a — + is taken and the "small A" approximation is 



used 5 , retaining only linear terms in A, Eq. (19) reduces to 



AP ~ 



-16A(l/2-|f/ M 3| 2 )|f/ e3 | 2 |f/ M 3| 2 (Asin2A-2sin 2 A) 



(20) 



in agreement with the result first reported in Ref. [25]. In this approximate limit, though not in 
general, AP^ is directly proportional to |?7 e3 | 2 (l/2 — |t/ M3 | 2 ). The more general dependence on 
s 13 and D is, of course, hidden in Eq. (19). 

The last-mentioned a — > limit may in fact be more relevant than being a mere device to 
get to Eq. (20). Given (1) the accidental feature that the atmospheric v^/ve ratio would have 
been ~ 2 in the relevant energy range if # u T oscillations were absent and (2) the fact of a 
near- maximal 623 , the effect of the solar neutrino squared mass difference Am 21 on atmospheric 
neutrino oscillations can be shown to be small [26] . It may therefore be useful to display also the 
analytical expression for the survival probabilities of muon neutrinos and antineutrinos with the 
assumptions of only the constancy of p and the neglect of O(a) terms. Taking the limit a — > 
and using the notation of Eqs. (16) and (18), one then has 



lim Ai 

31 -> 



1 



Am 2 31 (A+l) - (Am 



31/ 



5 From Table 1, the linear A approximation would appear to be accurate to within 90% only for neutrino energies 
below 3 GeV and pathlengths below 9000 km. 
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lim A2 = 0, 
lim A 3 = - 

a - 2 



2 \M 



Am^(i + l) + (A m ^) 



1/2 



for neutrinos with 

(Amj.f = AmlKA - 1)| [l + 4 S 2 3 (A - l)" 2 
and replacing A by —A for antineutrinos. Furthermore, it follows that 



lim 

a -> 


C/ M l 


— I s 23 s l3l 


lim 

a -> 


f/ M| 
<->2 1 


= c 23 , 


lim 

« -» 


U M \ 


— l S 12 c 13l 



(21) 

(22) 

(23) 
(24) 
(25) 



where sfl = sin^, cff = cos0;ff and 0^ is given in terms of 9 13 and A by the relation 



( A 



m 



31 



sin 2 20# = sin 2 20 i3 , 

13 {(Amir, 



sin 2 20 i3 



~ (A-cos20 13 ) 2 + sin 2 20 13 ' 
Eqs. (20) - (25) enable us to rewrite Eqs. (12) and (13) in the limit a — > as 

Q hm o P^(L) = 1 - P^(L) - P 2 ,(L) - P^(L), 



P^(L,A^-A), 



with 



Pi(L) = sm 2 « sin 2 20 23 sin 2 ^^^^ L, 



(26) 

(27) 
(28) 

(29) 



cos 2 0^ sin 2 20 23 sin 2 511 



P^(L) = sin 2 2^ sin 4 23 sin 



8P 



2 ( ^3/ 
4P 



(30) 



(31) 



We note that (L) vanishes in vacuum (when matter effects go to zero), but P 2 ^ 3 do not. In 
§3, we shall numerically study the behavior of P^ 2,3 (L) as functions of the neutrino energy E for 
different baseline lengths L. The expression for sin 2 20^f in Eq. (26) shows the effect of the MSW 
resonance when A equals cos 29 13, i.e. 



E = E r , 



Am| x cos 20i3 
2V2G F N e 



(32) 
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Though E vcs , of Eq. (32) has no explicit dependence on L, in practice an implicit L dependence 
creeps in through N e for path lengths involving significant spatial variations of the earth's density. 
Furthermore, at or near the resonance, Eq. (19) would not be trustworthy since, with s 13 small, A 
is rather close to unity - making some of the RHS terms blow up. At E — E xes ., sin 2 29^ reaches 
its maximum value, i.e. unity. For Am| x = (2 — 3) x 10~ 3 eV and a small #13, the resonance 
could occur [27] for atmospheric muon neutrinos passing through the earth with energies between 
5 and 10 GeV. At resonance, the conversion probability of a muon neutrino becomes sizable even 
for a small 9 13 . The effect of a small nonzero 9 13 could then be observed as an excess of upward 
going electron neutrinos. 

3 Numerical results on survival probability in matter 

We now present our results for the muon neutrino survival probability in vacuum and in 
matter. To obtain these, we have numerically solved the three-generation differential equation of 
motion of the neutrinos exactly, using the PREM profile for the earth matter densities [22] . Fig. 
1 shows PfA/j, as a function of L, for six different fixed values of E. This plot (as well as subsequent 
plots presented in this paper except where specified otherwise) has been generated with chosen 
benchmark values of the other concerned oscillation parameters Am^, Am^, sin 2 9 12 , sin 2 # 13 and 
5 - as tabulated in Table 2 - and two values of sin 2 6*23 as explained in the caption. The following 
features are evident from Fig. 1. 



Am 2 x = 2 x 10^ 3 


eV 2 


Am 2 21 = 8 x 10~ 5 


eV 2 


sin 2 12 = 0.28 




sin 2 20 13 = 0.1 




5 = 





Table 2: Chosen benchmark values of oscillation parameters, except 6*23, with assumed normal 
mass ordering. 



• Matter effects inside the mantle begin to be significant for E ^ 4 GeV, L ^ 4000 km and 
are largest for E between 5 and 6 GeV, reducing a bit for E = 7 GeV. 

• For E 4 GeV and inside the mantle, the survival probability in matter at the peaks 
invariably reduces from its value of near unity in vacuum. This effect increases as sin $23 is 
increased from 0.36 to 0.5 and beyond. 
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• For E ^ 4 GeV and inside the mantle, matter effects tend to increase the survival probability 
at the troughs from its vacuum value of = 1 — Ac\ z c\ z (l — c\ z c\ z ). This increase however 
reduces as the value of sin 2 #23 is raised from 0.36 to 0.5 and beyond. 

• For E ^ 4 GeV and L large but within the mantle (specifically between 8000 km and 10000 
km), the matter contribution to the survival probability, i.e. P MM — P^/f, can change sign. 

• The survival probability in matter has a fairly sharp discontinuity at the mantle-core bound- 
ary and behaves very differently with distance inside the core as compared with the mantle. 

For a given neutrino energy E, matter effects begin to be significant only when the neutrinos 
pass through matter densities which are close to their resonant value. Lower energy neutrinos need 
higher matter densities to hit the resonance condition and vice versa, as can be seen from Eq. 
(32). Since, for upward going atmospheric neutrinos, the average density of the earth increases 
as L increases, E res (av) of Eq. (32) decreases with L. In Fig. 2 we show the resonance energy 
E res (av), calculated from Eq. (32) with the benchmark (Table 2) values of sin 2 2#i3 and Am^, as a 
function of the distance L travelled by an atmospheric inside the earth. For each L, the average 
earth matter density has been calculated using the PREM density profile. For L < 1000 km, the 
resonance energy E res clearly exceeds 9 GeV. Since the flux of multi-GeV atmospheric neutrinos 
falls very fast [28] with energy and since we restrict our analysis to fully contained neutrino events 
in the detector with E < 11 GeV, we do not expect much matter effect for these short baselines. 

In order to understand the extent of matter effects for L > 1000 km, we need to take into 
account another factor apart from the difference between E and E res . That is the role of the 
Am^-dependent oscillatory terms in Eq. (27). In this context, it is instructive to start by 
looking at the extrema of the oscillatory term sin 2 (1.27Am 2 1 L/i?), where Am^ is in eV 2 , L is 
in km and E in GeV. We define -Espmini, -Espmax and -Espmin2 as the respective values of the 
energy corresponding to the first 6 minimum, the first maximum and the second minimum in the 
survival probability P w as L/E increases. In other words, these respectively correspond to the first 
maximum, the first minimum and the second maximum as L/E increases in the flavor conversion 
probability P^. Their values in vacuum are, 

P=e,r 

Pspmini = -(1.27)Arr&L, (33) 



Pspmax = -(L27)Am 2 1 L, (34) 

71 

£spmin2 = |-(1.27)Am 2 1 L, (35) 

are plotted, along with E res , against L in Fig. 2. We deduce from this figure that E res ~ -Espmini 
for L ~ 4500 km, while E res ~ -Espmax for L ~ 7400 km and E res ~ £'spmin2 for L ~ 9700 km. 



6 This means first or second (as the case may be) starting from the left in Fig. 1, but first or second starting 
from the right in Figs. 3 and 4. 
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Concentrating on the region 1000 km < L < 4000 km, we find that here 7.4 GeV < E res < 9.6 
GeV, and E res ^> E S pmini so that once again the survival probability P M(U is not much affected 
by earth matter. In contrast, for L > 4000 km and with E between 5 and 7 GeV, large matter 
effects ensue, as evident from Fig. 1. 

The three terms P* P 2 M and P 3 as defined in Eq. (27) and given analytically in Eqs. 
(29)-(31), are shown in Fig. 3 as functions of the neutrino energy E for the three special baseline 
lengths mentioned above: L = 4500 km (upper panels), 7000 km (middle panels) and 9700 km 
(lower panels). The dotted blue lines in Fig. 3 show these three terms in vacuum for sin 2 6> 2 3 = 1/2, 
while the other lines show them in the presence of matter for different values of sin 2 8 2 s, as explained 
in the caption. For these plots we have taken Am| x = 0, while the other parameters are as in Fig. 
1. We discuss below the behavior of the three terms P^ 2,3 separately 7 . 

1. P* is proportional to sin 2 6^, cf. Eq. (29). Though, approximately equal to sin 2 #13 for 
E <C E res , sin 2 9fl increases with E and reaches the value 1/2 at E = E res . Beyond E res , 
sin 2 9^ continues to rise with E till it saturates to the maximum value of unity for E S> E res . 
For L = 4500 km, one has 5_ = sin 2 [{Am^(A + 1) - (Am^) M }L/8P] ~ 1 for most part of 
the energy range in Fig. 3. However, for L = 7000 km, S- is ~ 0.2 in magnitude at E — 1 
GeV and rises monotonically to ~ 0.9 at E = 9 GeV. As a result, for both baseline lengths 
of 4500 km and 7000 km, P^ increases monotonically with energy. For L = 9700 km, S- is 
~ 0.8 around E = 1 GeV, but falls to almost zero at E = 7 GeV and rises yet again to ~ 0.3 
at E = 9 GeV. The consequence is that P* which is a product of a continuously rising 
sin 2 9™ and the oscillatory S- factor, has a complicated behavior as a function of energy at 
L = 9700 km. We recall here that = AP^ = P^ - P^ ac is a pure matter effect term 
and hence always positive. 

2. The second term P 2 ^ is the dominant term in vacuum. For E <C E res , there is very little 
matter effect in it and P 2 ^ remains close to its vacuum value. As E increases towards E res 
and beyond, cos 2 #{*f, which appears as a factor in P 2 ^, decreases towards 0.5 at E res and 
eventually to near-zero for E S> E res . For E < E res and all values of L, the matter-dependent 
oscillatory factor S + = sin 2 [{Am| 1 (A + 1) + (Am 2 1 ) M }L/8P] closely follows the behavior of 
the corresponding sin 2 (Am 2 1 L/4P) term in vacuum. But the situation changes drastically 
for E > E res . Now the said S + factor is almost unity at L — 4500 km and 9700 km, but 
nearly vanishes at L — 7000 km. Thus we can describe all the cases for P 2 ^ in the following 
way. While P 2 ^ follows its vacuum oscillation pattern for E < E res , it falls to almost zero 
as E increases beyond E res . The latter behavior ensues at L — 4500 km and 9700 km due 
to the fall of cos 2 6>ff to zero, whereas at L — 7000 km it is caused by S + nearly vanishing. 
For E ~ E res , AP 2 ^ = P 2 ^ - P^ ac is positive at L = 7000 km but negative at L = 4500 
and 9700 km. 

3. Turning to P 3 ^, we see from Eq. (31) that its proportionality to sin 2 20^ sin 2 [(Am 2 1 ) M L/4P] 
makes it very small in vacuum owing to the smallness of sin 2 20 13 . In matter, however, the 
growth of sin 2 29f^ with energy (for E ^ E res ) till its maximum value at E — E res makes 

7 For a similar discussion on the dependence of matter effects in P w on L and E see [21]. 
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P 3 M increase appreciably in this energy range. For our example cases shown in Fig. 3, 
(E ~ 5 GeV and maximal # 23 ), it increases from about 0.01, 0.0 and 0.025 in vacuum at 
L = 4500, 7000 and 9700 km respectively to about 0.14, 0.20 and 0.25 in earth matter. This 
means that, at E ~ E res , AP 3 ^ = P 3 ^ — P 3 ^ ac is in general positive for all L. 

The impact of earth matter effects is largest for any L at E = E res when the sum of the three 
terms P* P 2 ^ and P 3 M becomes most different from its value in vacuum. In fact, P^ ac (P^ ac ) is 
exactly (nearly) zero for all E, while - for E ~ 5 - 7 GeV - Pjjfi™ yields SPMIN1 at L = 4500 km, 
SPMAX at L = 7000 km and SPMIN2 at L = 9700 km. While matter effects raise P 1 ; 3 for all L 
and E (cf. Fig. 3), apart from mildly increasing P; M at SPMAX for L = 7000 km, they decrease 
it when L = 4500 (SPMIN2) and 9700 km (SPMIN2). There is thus an important difference 
between the SPMAX and the two SPMIN cases: matter effects come with the same relative sign 
for all the three terms at SPMAX, but increase P^ 3 and decrease P^ at SPMIN1 and SPMIN2. 
The SPMAX case at L — 7000 km thus exhibits the largest effects of earth matter in the total 
survival probability P M(U for E ~ 5 — 7 GeV. In contrast, these effects partially cancel in the two 
SPMIN cases between the increase in P^ 3 and the decrease in P% For L = 4500 km (SPMIN1), 
this cancellation is particularly important and the total P MM has very little residual matter effect. 
Therefore we shall ignore SPMIN1 in our subsequent matter effect discussion. For L = 9700 km 
(SPMIN2), the cancellation is not complete and appreciable matter effects do persist in P MM . As 
noted earlier, one characteristic of the SPMAX case (L = 7000 km and E ~ 5 GeV) is the decrease 
in Pjxjx from its maximal value of near-unity in the presence of matter since all the three terms 
P^ 2 ' 3 increase. In contrast, P MM increases on account of matter in the SPMIN2 case (L = 9700 
km, E ~ 5 GeV, as exemplified here) because of the dominant decrease in P 2 ^. 

We can comment on the respective roles of 13 and # 2 3 in determining the effects of matter 
on the muon neutrino survival probability. First, the value of #13, of course, directly controls 
the extent of such effects. Figs. 1-3 have been generated with sin 2 26>i 3 = 0.1. For the same 
L (and hence the same average matter density), the value of E res is larger for a smaller #13, 
cf. Eq. (32). On the other hand, the value of (Amj^ is smaller in that case for a constant 
matter density independent of L. The result is a larger mismatch between E res (at which sin 2 29f^ 
becomes maximal) and the energy where P 2 ^ 3 reach their extremal values. Matter effects in P MM 
are therefore reduced for a smaller 9 13 . Second, Fig. 1 carries evidence of the dependence of the 
contribution of matter effects to P w on sin 2 #23 . The nature and extent of this dependence is 
more clearly brought out by Eqs. (29-31) and Fig. 3. While P^ 2 are seen to be proportional to 
sin 2 2#23, the third term P 3 ^ goes as sin 4 623- A noteworthy point in this connection is that, when 
623 is close to being maximal, any change in its value affects sin 2 2# 2 3 much less than sin 2 6*23. For 
instance, when #23 decreases from 45° to 40°, sin 2 2# 23 changes from 1 to 0.97, i.e. by 3%, while 
sin 2 # 23 is reduced by 16% from 0.5 to 0.41. This is reflected in essentially no visible change in 
P^ 2 when sin 2 #23 is reduced 8 (increased) from 0.5 to 0.4 (0.6), cf. Fig. I. On the other hand, 
P 3 does increase appreciably as sin 2 #23 is increased from 0.4 to 0.5 leading to an increased #23- 
sensitivity in a muon neutrino disappearance experiment on account of matter effects. As noted 
earlier, the net decrease (increase) of the total survival probability P MM due to matter effects at 
SPMAX (SPMIN2) increases (decreases) as the value of sin 2 #23 is increased. This behavior can 

8 sin 2 2#23 is the same for sin 2 #23 = 0.4 or 0.6. 
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now be understood from the dependence of P^ on sin 2 #23 ■ Let us first consider SPMAX. Since 
the increase due to matter effects in P^ causes a decrease in the total survival probability P MM 
at SPMAX and since in matter increases as sin 2 #23 increases, the matter-induced decrease of 
P MM at SPMAX increases with increasing sin 2 823 ■ For the SPMIN2 case, the dependence of on 
sin 2 023, as discussed above, is more subtle. In reality, the main role of the increase in P£ due to 
matter is to wash out the corresponding decrease in P 2 ^. Therefore, the matter-induced relative 
increase of the survival probability P MM at SPMIN2 decreases as sin 623 increases. 

It is pertinent to stress a couple of points here. First, the strong dependence on sin 2 6> 23 of the 
total survival probability P w in matter coming from P 3 , as discussed above, implies that observed 
atmospheric v^, events could be used not only to probe \D\ = 1 1/2 — sin 2 #23] but also the sign 
of D to determine if #23 < 7r/4 or > 7r/4, thus resolving the ambiguity regarding the octant of 
#23- The second issue is the effect of earth matter on the shape of the survival probability as a 
function of the neutrino energy. This is shown in Fig. 4 for six different choices for L. The various 
choices of sin 2 #23 are given in the figure caption while the other parameters are the same as in 
Fig. 1. This figure again tells us that matter effects lower (raise) P M/i from its vacuum values at 
SPMAX (SPMIN1 and SPMIN2), resulting in a relative change of the shape of P MM as a function 
of E. These matter-dependent effects will of course be projected onto the observed atmospheric 
muon neutrino spectrum which will come folded with P^. A detector with good energy resolution 
should be sensitive to the spectral shape of atmospheric muon neutrinos and hence should pick 
out the shape distortion owing to matter effects. 

4 Up-down asymmetry in an ICAL-like detector 

A very powerful variable, which clearly displays oscillation effects for atmospheric neutrinos (an- 
tineutrinos) , is the up-down ratio [29] U N /D N (Ua/D a ). Here U N (C/4) is the number of events 
recorded for the "upward" muon neutrinos (antineutrinos) coming with cos£ < and D N (Da) 
is the number of atmospheric events recoded for the "downward" muon neutrinos (antineutrinos) 
with cos£ > 0, £ denoting the zenith angle of the v^jv^ trajectory. Downward (upward) neutrinos 
and antineutrinos cover short (long) distances. Hence downward ones hardly undergo any flavor 
transformation while upward ones are subjected to full flavor oscillations, the relevant oscillation 
lengths being in the range 10 2 — 10 3 km. Thus the up-down ratio U/D yields crucial information 
on the degree of oscillations for upward neutrinos and antineutrinos. Being relatively insensitive 9 
to the uncertainties in the atmospheric neutrino flux calculations, it is a reliable measure of the 
survival probability P M/i . As such, it can be used to study both the vacuum oscillations of muon 
neutrinos and the effect of matter on them [20]. We shall consider the above ratio for atmospheric 
muon neutrino and antineutrino events in a large magnetized iron calorimetric detector. Such a 
detector named ICAL, has been conceived in the proposed underground neutrino laboratory INO 
[18] in India. While the initial proposal is for a detector mass of 50 kton for ICAL, the aim is 
to enlarge it to a final detector size of 100 kton. The proposed detector would have a modular 
structure with stacks of ~ 6 cm thick magnetized iron plates interleaved with ~ 2.5 cm of resistive 
plate chambers (RPC) made of glass as the active detector material. A uniform magnetized field 

9 We have more to say on this issue later in this section. 
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of 1 — 1.4 Tesla would enable the charge discrimination of muons distinguishing between v^- and 
v^— induced events. The full detector would be divided into three modules, the modular structure 
allowing the start of detector operations with those modules that are ready even when others are 
still under construction. Further details are given in [18]. 

The expected number of /i (or p) events, induced by oscillating atmospheric muon neutrinos 
(or antineutrinos) and recorded in a detector such as ICAL, is given by 

Ni= n T / dcosM dgRfag) dE dE'R(E,E') a^E') ei(E') 

(36) 



dE' dcos? " ( ' * } + dE'd cos? Pl ' l{E ' * } 



where, E' and E are the true and experimentally reconstructed neutrino energies respectively 
and R(E, E') is the energy resolution function of the detector. Likewise, £' and £ are the true 
and reconstructed neutrino zenith angle and -R(£, £') is the zenith angle resolution function of the 
detector. It is convenient for us to adopt 10 Gaussian forms for the energy and length (zenith angle) 
resolution functions of the detector: 

^-^kM^f) • (37) 

where x is E or £ (L). For the energy resolution we assume 2a e> = 0.3E'. The detector resolution 
for the zenith angle and hence the distance travelled by the neutrino is expected to be somewhat 
better and we assume 2a L > = 0.2L'. These numbers are more or less the same as the full widths 
of the nonGaussian resolution functions given in [18]. The total distance L and the zenith angle 
£ are related by 



(R e + h) 2 - R e 2 sin 2 £ - R e cos£ , (38) 

where R e is the radius of the earth and h the height of the atmosphere. Among the rest of the 
quantities of Eq. (36), tit denotes the number of target nucleons times the total live time of the 
detector, ai(E') is the total viN — > IX scattering cross-section of a v\ of energy E' and ei(E') is the 
trigger efficiency of the magnetized calorimeter. Here the index / can be \x or Jx and d 2 <pi/dE'd cos £' 
is the differential flux of atmospheric neutrinos v x . Finally, Pu gives the survival probability for 
the atmospheric (or v^) and Pyx the transition probability of v\ to u h V being either e or e 
correlated to / being li or /L 

We shall present our results for a 100 kT of active detector mass and ten years of running, as- 
suming a conservative 50% trigger efficiency 11 for the detector 12 . In order to calculate the number 
of detected events, we use the latest three-dimensional atmospheric neutrino fluxes provided by 
Honda et. al. [28], which is also used by the INO collaboration [18]. For the reaction cross-section 

10 We have checked that for the E and L bin sizes used in this paper the resolution functions make only little 
difference to the final results. Hence, the precise form of the resolution function is unimportant. 
11 In reality, this could be as high as 80%, S. R. Dugad, private communication. 
12 For any other detector mass and efficiency, the statistics can be accordingly scaled. 
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we consider only the DIS process and use the cross-sections given by the CTEQ collaboration [30] . 
We then distribute the number of fj,~ events induced by the (P M ), into various zenith angle 
and energy bins. In the range E — 1 — 11 GeV, the events are divided into five energy bins, each of 
width 2 GeV. The zenith angle binning is presented in Table 3. The bins 1 — 6 (7—12) correspond 
to cos£ > (< 0) and hence contain downward (upward) going neutrinos/antineutrinos. Upward 
going ones have been binned according to L m , the distance they travel within earth matter. The 
range of L m (in km), corresponding to each zenith angle bin, is given in Table 3. The bins for 
downward neutrinos are arranged to have a structure similar to that of upward neutrinos. We 
simulate the "data" for the prototype detector, assuming certain "true" values for the oscillation 
parameters namely those benchmarked in Table 2, and distribute them into the five energy bins 
and the twelve zenith angle bins. Thus we have 5 x 12 = 60 bins of data for the muon events. 
The iron detector we consider is magnetized and therefore has charge discrimination capability. 
So we have 60 bins of /i~ data for the neutrino channel and another 60 bins of fj, + data for the 
antineutrino channel. Our total data set therefore comprises 120 bins. 

The up-down ratio Un/Dn for z/ M 's is shown in five different energy bins in Fig. 5. However, 
instead of integrating over all upward and downward zenith angles, we have chosen to divide 
the said ratio into six zenith angle bins. Each of the panels in Fig. 5 contains Un/Dn for a 
certain range of the modulus of the zenith angle as shown in the legend. For instance, the first 
(last) panel on the top left (bottom right) contains the ratio Un/Dn for < |cos£| < 0.157 
(0.785 < | cos£| < 1) and calculated by taking upward going neutrinos with —0.157 < cos£ < 
(— 1 < cos£ < —0.785) and downward going neutrinos with < cos£ < 0.157 (0.785 < cos£ < 1). 
The black and magenta solid lines in Fig. 5 describe the Un/D n spectrum for the realistic case 
of upward neutrinos travelling through earth matter and with sin 2 6*23 = 0.5 and sin 2 #23 = 0.36 
respectively. In contrast, the similar dashed lines show the hypothetical Un/Dn spectrum for 
the corresponding values of sin 2 #23, if the upward neutrinos were assumed to be travelling in 
vacuum even inside the earth. For the same value of sin 2 #23, a comparison of the solid line with 
the corresponding dashed line shows the effect of earth matter in changing the muon neutrino 
up-down ratio. Furthermore, the degree of this change is seen to depend on the value of sin 2 #23 . 

An even better way of bringing out matter effects 13 in atmospheric muon neutrinos is to look 
at the difference Un/Dn — Ua/Da, where the subscript A now refers to P M 's. The up-down ratio 
U/D yields P^ which is the same for neutrinos and antineutrinos in vacuum so long as CPT is 
conserved. Therefore, the difference in the up-down ratio between neutrinos and antineutrinos 
gives a direct measure of the matter effects in P^. In fact, it relates to AP W of Eq. (19). Fig. 
6 shows our results for the same choice of zenith angle and energy bins and for the same set of 
oscillation parameters with an assumed normal mass ordering as made in Fig. 5. Substantial 
matter effects in terms of the deviation of the difference Un/Dn — Ua/Da from zero are con- 
cretely shown in Fig. 6 and in a more pronounced manner in the right hand panels. Once again, 
the strength of the matter part in the survival probability is seen to depend on the value of sin 2 #23 . 

We can summarize the main characteric features of Figs. 5 and 6 as follows: 

13 For a normal (inverted) mass ordering, matter effects at baseline lengths L > 1000 km for v^s (P^'s) are 
significantly larger than for P M 's (*v's), cf. discussion at the end of the section. 
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• For a normal mass ordering (Am^ > 0), atmospheric muon neutrinos undergo significant 
matter effects for E ^ 3 GeV and L ^ 4000 km. A similar statement holds true for muon 
antineutrinos if Am^ < 0. 

• Matter effects increase (decrease) the ratio Ujy/D N in the E and cos£ bins corresponding to 
the SPMINs (SPMAX) for a normal mass ordering. For an inverted mass ordering, a similar 
statement can be made about Ua/Da- 

• The largest impact of earth matter comes in the bin E = 5 — 7 GeV and —0.628 < cos£ < 
-0.471 (6000 < L m < 8000 km), corresponding to the SPMAX. Here U N /D N - U A /D A 
reduces for a normal mass ordering to —0.3 to —0.25, depending on the value of sin 2 #23. 

• The dependence of the strength of matter effects on the value of sin 2 #23 is most clearly 
brought out in the energy and zenith angle bin mentioned in the preceding bullet. Reducing 
sin 2 #23 from the maximal 0.5 to 0.36 brings a nearly 10% change in the difference Ujy / 'Djy — 
Ua/Da- This highlights the possible role of matter effects in achieving a better sensitivity 
to sin 2 823 and hence a better resolution of its difference from maximality. 

Let us make two further comments. First, the main advantage of the up-down asymmetry 
parameter, constructed from directly recorded events, is its insensitivity to the error in the absolute 
normalization of the atmospheric neutrino flux, on account of cancellations between the upward 
and downward fluxes. Recall that this uncertainty has been the main source of systematic error in 
extracting the mass squared difference and mixing parameters from any set of observed atmospheric 
neutrino events. Additionally, the charge discrimination capability of an ICAL type of a detector 
should be able to bring out earth matter effects in the atmospheric neutrino signal much more 
effectively by constructing the difference U^/D^ — Ua/Da from recorded events. The second 
comment relates to our assumption of a normal neutrino mass ordering, automatically implying 
larger matter effects in the rather than the channel for L > 1000 km. In case Nature has 
chosen an inverted ordering of the neutrino masses, the above features would hold qualitatively 
except that the situation would be reversed between the and v^. In particular, the major change 
in the plot of Fig. 6 would be that of a sign. 



Our twin goal is to check the sensitivity of atmospheric v^s and P M 's to any deviation of 6*23 
from its maximal value 7r/4 and to see the effects of earth matter on this sensitivity. In order to 
achieve this end, we perform a statistical analysis of simulated "data" generated in an ICAL-like 
calorimeter, assuming certain "true" values for the oscillation parameters. For our error analysis 
we define a x 2 parameter following [31, 32, 33] as 



5 Details of the Statistical Analysis 




2 



(39) 
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bin 


cos£ 


Total distance L in km 


Distance in earth L m 


1 


1.000 to 0.785 


15.0 


- 19.1 





2 


0.785 to 0.628 


19.1 


- 23.8 





3 


0.628 to 0.471 


23.8 


-31.7 





4 


0.471 to 0.314 


31.7 


-47.3 





5 


0.314 to 0.157 


47.3 


-91.5 





6 


0.157 to 0.000 


91.5 - 


- 437.4 





7 


0.000 to -0.157 


437.4 - 


- 2091.5 


- 2000 


8 


-0.157 to -0.314 


2091.5 


- 4047.3 


2000 - 4000 


9 


-0.314 to -0.471 


4047.3 


- 6031.7 


4000 - 6000 


10 


-0.471 to -0.628 


6031.7 


- 8023.8 


6000 - 8000 


11 


-0.628 to -0.785 


8023.8 - 


- 10019.1 


8000 - 10000 


12 


-0.785 to -1.000 


10019.1 


- 12757.0 


10000 - 12500 



Table 3: The zenith angles bins used. 



In Eq. (39) N^ ata are the observed number of events in bin n and a^ at are the statistical errors. 
The systematic errors in the data and the theory are accounted for through the set of pulls 
The latter are defined in such a way that the number of expected events N^ heory , in bin 
n, corresponds to = and their la deviations are given by = ±1. The theoretical and 
experimental uncertainties then cause the expected number of events to shift to jV* heor?/ : 



jytheory j^theory 



+0(&. (10) 



k=l 



In Eq. (40) J /V* /ieor?/ has been expanded in powers of keeping only linear terms. The quantities 
give the fractional rate of change of N^ heory due to the fcth systematic uncertainty. 
The most important theoretical systematic uncertainties come from our lack of knowledge of 
the predicted atmospheric neutrino fluxes. In order to cover these, we take an absolute normal- 
ization error of 20%, a "tilt" factor of 5% to account for the spectral uncertainty, an uncertainty 
of 5% in the neutrino-antineutrino flux ratio and an uncertainty of 5% in the zenith angle de- 
pendence. For the DIS cross-sections, we consider an uncertainty of 10% - assuming a modest 
improvement in our present understanding of these cross-sections by the time ICAL starts oper- 
ating. For the other experimental systematic uncertainties, we take a 5% consolidated error. We 
take these somewhat arbitrarily assumed plausible systematic errors for the proposed experiment, 
pending actual estimates to be provided by the INO collaboration in future. We calculate the 
uncertainties in Eq. (40) using an approach similar to that used in [31]. Therefore, ir* is taken 
as 20%, 10% and 5% of the number of events for all bins when k corresponds to the absolute 
flux normalisation, cross-section uncertainty and experimental systematic error respectively. The 
7r^s for k corresponding to the "tilt" factor and zenith angle dependence of the and flux 
are calculated using the method detailed in [31]. For the flux uncertainty coming from neutrino- 
antineutrino flux ratio, we take 7r^ = 2.5% of the number of events for all the neutrino bins and 
7r^ = —2.5% of number of events for all the antineutrino bins. 
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Unless otherwise stated, for the numerical analysis done in this paper, the number of "data" 
events N% ata in the nth bin is simulated at "true" values 14 of the parameter set {p (true)}, p 
covering Am^, Am^, sin 2 9 12 and 5, which are set at the corresponding benchmark values listed 
in Table 2. The choice of a normal mass ordering is for the sake of being definite rather than 
on account of any theoretical prejudice. We shall also present some results for the inverted mass 
ordering. The true values of the remaining oscillation parameters sin 2 #23 and sin 2 #13 will be varied 
to account for their impact on matter effects and will be stated explicitly as and when used. (Any 
change, whenever made in the assumed true values of the other parameters, will also be mentioned 
explicitly.) For each such simulated data set, we do a statistical analysis using Eqs. (39) and (40) 
to find the sensitivity of ICAL to any deviation of #23 from maximality as well as to the fixation of 
its right octant. A^ /ieor?/ is calculated for a given set of oscillation parameters using Eq. (40). The 
RHS of Eq. (39) is then minimized with respect to the pulls to obtain xltm as a function of the 
oscillation parameters. The latter function is then further minimized by varying the oscillation 
parameters within their allowed ranges to obtain the sensitivity plots for 823 ■ 

The extent of earth matter effects is determined by the true value of sin 2 #13. Therefore, we 
also include the constraints on sin 2 9 13 expected to ensue from the combined data expected from 
the next generation reactor [34] and long baseline accelerator [35] experiments. So we define the 
combined x 2 for the ICAL experiment as 



In Eq. (41), xltm * s as gi ven by Eq. (39). Furthermore, the second RHS term tries to take 
into account the above-mentioned bounds on sin 2 # 13 from future laboratory experiments with the 
denominator a s 2 denoting the la uncertainty in it which is assumed [36] to be at the level of 



6 Deviation of #23 from its maximal value 

Let us turn now to the focal point of this paper: how well will the deviation of #23, if any, from 
its maximal value be probed with atmospheric neutrinos in the foreseeable future? We take the 
combined information from 1 MtonY of simulated data in ICAL, adding on the constraint on 
sin 2 9 1 3 to come from future reactor and accelerator experiments. Our procedure is to generate 
the data at a certain nonmaximal value of sin 2 6*23 (true) and then fit this data with the maximal 
muon neutrino mixing angle 6*23 = 7r/4, choosing different values of sin 2 6*23 (true). The results are 
displayed in Fig. 7. At each point in this Am 2 1 (true) - sin 2 $23 (true) plane, we simulate the 120 
bin data in ICAL, taking sin 2 #13 (true) = 0.00 (left panel), 0.02 (middle panel) and 0.04 (right 
panel) respectively and assuming for all the other parameters the benchmark values given in Table 

14 In what follows, wc shall always distinguish between the "true" values of the oscillation parameters chosen by 
Nature at which we generate our projected data and the fitted values of those parameters which are constrained 
by the same data. The true values will henceforth be written as sin 2 #23 (true) etc., while the fitted values will be 
referred to as sin 2 #23 etc. 




(41) 



3.5% 
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2 to be true. Each such data set is then fitted back using Eqs. (40), (39) and (41) with N l ^ eorv 
calculated for maximal #23 mixing. The parameters Am^ and sin 2 #13 are allowed to take any 
possible values in the fitted j\T^ eon/ , while Am| x and sin 2 12 are allowed to vary freely within 7% 
and 15% of their assumed true values respectively 15 . 

The upper and lower panels of Fig. 7 correspond to respective exposures of 1 MtonY and 3.37 
MtonY 16 . The regions of sin 2 #23 (true) and Am| 1 (true) within the white, blue and green bands 
of Fig. 7 show the true values of those quantities for which the distinction of a maximal from a 
true nonmaximal value of #23 will not be possible at the la, 2a and 3a levels respectively for the 
specified exposure. The broken lines give the corresponding limits of sin 2 #23 (true) in case earth 
matter effects were deliberately switched off by hand. Among those the long-dashed, dot-dashed 
and dotted lines respectively yield the la, 2a and 3a limiting values of sin 2 #23 (true). A comparison 
of the broken lines with the corresponding continuous lines show the following feature. Matter 
effects tend to increase somewhat the sensitivity of ICAL to test the maximality of sin 2 #23 • 

We present in Table 4 the intervals of D(true) = 1/2 — sin 2 #23 (true) beyond which max- 
imal mixing could be ruled out at the 3a level, with and without matter effects. The latter 
corresponds to a fictitious environment of pure vacuum and has been included to see the quanti- 
tative role of matter effects here. Specifically, for Amf^true) = 2.0 x 10~ 3 eV 2 and sin 2 #13 (true) 
= 0.04 (0.00), sin 2 9 23 (true) can be distinguished by ICAL from the maximal value of 0.5 at the 
3a level within 17% (20%) from 1 MtonY of exposure and within 11% (14%) if the statistics was 
increased to 3.37 MtonY. The corresponding ranges of D(true) are what Table 4 lists. This is 
comparable to the sensitivity of the combined data from the forthcoming accelerator-based long 
baseline experiments 17 to a deviation from maximality of sin 2 #23, which is [12] ~ 14% at 3a (for 
Am| 1 (true)= 2.5 x 10~ 3 eV 2 ) 18 . Our sensitivity to D is also comparable to that expected with 
atmospheric neutrinos in very large futuristic water Cerenkov detectors. For statistics that is 20 
(50) times the current SK statistics, denoted as SK20 (SK50), a very large water Cerenkov atmo- 
spheric neutrino experiment is expected to test a deviation from a maximal sin 2 #23 upto [9] 23% 
(19%) at 3a. 

Let us make some observations here. While atmospheric neutrino data in a very large water 
Cerenkov detector could test [9] maximality in sin 2 #23 almost independently of the values of 
sin 2 #13 (true) and Am 2 1 (true), the sensitivity at a detector like ICAL to D would depend on both 
these oscillation parameters. The dependence of the latter on sin 2 #13 (true) comes from the greater 
importance of matter effects here. The impact of Amf^true) on the measurement D comes from 
the dependence of sin 2 #23 sensitivity to the spectral shape of the data at ICAL, cf. §3. For a 
larger Am| 1 (true), there is a larger averaging of the oscillation signal and the sin 2 #23 sensitivity 
is reduced slightly. Such a fact is evident from Figs. 8 and 9 which display the muon (/i~) zenith 

15 These are the bounds on the solar parameters expected from the future solar and long baseline reactor neutrino 
experiments. See for example [37] for a recent detailed discussion. 

16 3.37 MtonY would roughly be the statistics needed in ICAL to match the number of fully contained muon 
events in the experiment SK20 which is described later in this section. 

17 The authors of [12] use a combination of simulated data set from five years of running of MINOS, ICARUS, 
OPERA, T2K and NO^A each, expected to come in the next ten years. 

18 This estimated uncertainty will be quite different if the true value of Ato^ deviates substantially from 2.5 x 10~ 3 
eV 2 . 
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Range of D(true) in units of 10 2 




matter 


vacuum 


sin 2 13 (true) 


0.00 0.02 0.04 


0.00 0.02 0.04 


1 MtonY 


[9.8,-9.7] [8.8,-10.3] [8.4,-10.7] 


[9.4,-9.4] [9.2,-10.2] [9.0,-11.1] 


3.367 MtonY 


[7.0,-7.0] [5.9,-7.2] [5.5,-7.2] 


[6.6,-6.7] [6.3,-7.3] [6.0,-8.0] 



Table 4: Simulated ranges of deviation D from maximal sin 2 # 2 3, allowed at 3a, after two different 
exposure times at ICAL. Here Am| 1 (true) has been taken to be 2.0 x 10" 3 eV 2 . 

angle spectrum for Am 2 x = 2.0 x 10~ 3 eV 2 and 4.0 x 10~ 3 eV 2 respectively. The four panels in both 
figures show the zenith angle dependence of the muon events in four different energy bins which 
are labeled. While the solid black and the magenta dotted lines are for sin 2 #23 = 0.5 (maximal) 
and 0.4 respectively, the blue dot-dashed lines correspond to unoscillated neutrinos. A comparison 
between Figs. 8 and 9 reveals an interesting feature. The case with the higher value of Am| x has 
a spectral distortion which is less than that with the lower value, especially when E is between 
3 and 5 GeV and when E is between 5 and 7 GeV, which are statistically the most important 
energy bins. The reason is that the oscillations are faster for Am^ = 4.0 x 10~ 3 eV 2 , leading to 
partial averaging and a smaller spectral distortion in the resultant signal. In fact, for the same 
reason, matter effects also are slightly more important for Am| x = 2.0 x 10~ 3 eV 2 , as compared 
to the higher value quoted above. Such an increased averaging, especially in the statistically more 
relevant lower energy bins, results in a slight fall in the sensitivity of ICAL to the precision of the 
oscillation parameters. This happens since the use of spectral distortion is a crucial component of 
the latter in relation to oscillation parameters, including the deviation D from the maximality of 
sin 2 6» 23 . 

7 Sensitivity to the octant of # 2 3 

We now come to the determination of the sign of D = 1/2 — sin 2 9 23 which decides whether 
O23 < ir / 4 (D positive) or 9 2 s > n / 4 (D negative), i.e. whether 623 lies in the first or second octant 
of the (0, 2tt] range. In a two-generation picture, the dependence of the survival probability P w 
only on sin 2 2# 2 3 makes it impossible to fix the octant of #23- Even in a three-generation framework, 
given that #i 3 and a are small, the dominant part of in vacuum is still found to depend only on 
sin 2 2#2 3 . Therefore, once again the fixation of the octant of #23 would be nearly impossible from 
a study of vacuum oscillations in surviving muon neutrinos. This degeneracy, however, is broken 
in matter owing to an additional strong sin 4 #23 dependence in P MM , as clear from the third term 
P 3 M in Eq. (31). Since we have already demonstrated the significant potentiality of the ICAL-like 
detector in deciphering the effects of earth matter on its observed muon neutrino/antineutrino 
events, we expect such a detector to be sensitive to the octant of 6*23. That is what we now 
numerically investigate. 

The simulated event spectrum at this detector is first generated at some nonmaximal "true" 
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value of 623. Then we try and see if the wrong octant of 6*23 can be ruled out at a good C.L. 
Here we use anticipated constraints on sin 2 6*13 from future reactor and long baseline experiments, 
but we have allowed all other parameters except sin 2 6*23 to vary freely in the fit. We consider 
three different cases of sin 2 #23 (true) values, namely 0.42, 0.46 and 0.54. The "true" values of 
the other oscillation parameters are as listed in Table 2 with an assumed normal mass ordering. 
Fig. 10 shows the results of our statistical analysis based on simulated data from 1 MtonY of 
ICAL exposure with the left, middle and right panels corresponding to the values 0.42, 0.46 and 
0.54 respectively of sin 2 #23 (true). In each panel the long dashed magenta, short dashed blue, 
dot-dashed green, dotted red and solid black lines show the y 2 for sin 2 # 13 (true) = 0.00, 0.01, 0.02, 
0.03 and 0.04 respectively. For every nonmaximal sin 2 6*23 (true), there exists a sin 2 #23 (false) which 
is given by 

sin 2 #23 (false) = 1 - sin 2 9 23 (true) (42) 

on the other side of ix/A. For a vanishing sin 2 #13 (true) there are no matter effects and the x 2 
corresponding to both the true and false values of sin 2 #23 are the same. Hence they are allowed 
at the same C.L. and one fails to fix the octant of #23 m this case. However, for sin 2 # 13 (true) 7^ 0, 
matter effects bring in an octant sensitivity and a false sin 2 #23 solution can be ruled out, provided 
D(true) is not too close to zero. For a given sin 2 #13 (true), the C.L. at which this can be done for 
our illustrative cases of sin 2 #23 (true) can be read out using Fig. 10, from the difference in the x 2 
between the true and false solutions. 

In order to obtain the limiting value of sin 2 623 (true) which could still allow for the determina- 
tion of sgn(D) we define 

Ax 2 = X 2 ( sm2 #23 (true), sin 2 # 13 (true), others(true)) 

— X 2 (sin 2 623 (false), sin 2 #13, others), (43) 

with 'others' comprising Am 2 l5 Am^, sin 2 9u and 5. These, along with sin 2 613, are allowed to 
vary freely in the fit. Eq. (43) gives us a measure of the C.L. at which sin 2 6*23 (false) is disfavored 
for a given sin 2 623 (true). Fig. 11 shows A% 2 as a function of sin 2 6* 2 3(true) for four values 0.01, 
0.02, 0.03 and 0.04 of sin 2 # 13 (true), corresponding respectively to the blue dashed, and green 
dot-dashed, the red dotted and the black solid lines. The left hand panel shows the results for a 
normal mass ordering of neutrinos while the right hand panel corresponds to the situation if the 
mass ordering for the neutrinos were inverted. The range of sin 2 6*23 (true), for which sin 2 #23 (false) 
can be ruled out at the 3a level, is visible from the figure. In particular, for a normal neutrino 
mass ordering, sin 2 9 2 3 (false) should be excludable at the 3a level from 1 MtonY of ICAL exposure 



for the following cases: 

sin 2 # 23 (true) < 0.361 or > 0.633 for sin 2 ^3 (true) = 0.01, (44) 

sin 2 # 23 (true) < 0.402 or > 0.592 for sin 2 ^3 (true) = 0.02, (45) 

sin 2 #2 3 (true) < 0.415 or > 0.580 for sin 2 0i 3 (true) = 0.03, (46) 

sin 2 e 23 (true) < 0.421 or > 0.573 for sin 2 9 13 (true) = 0.04. (47) 
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For an inverted mass ordering, the sensitivity of ICAL to sgn(D) is seen to be slightly less 
from Fig. 11. This is the result of the following fact. When the mass ordering is inverted, 
larger matter effects appear for u^s instead of i^'s at baseline lengths L > 1000 km. Owing 
to smaller interaction cross-sections for the antineutrinos, we roughly expect only half as many 
z/^-induced antimuon events in the detector as compared with i^-induced muon events. Therefore, 
the statistical power of the experiment for deciphering sgn(D) goes down by roughly a factor of 
half in this case. 

In order to quantify the importance of distinguishing the //^-induced muon events from the 
//^-induced antimuon events by means of the magnetic field in the detector, we have repeated our 
X 2 analysis by "switching off" the magnetic field - or in other words by taking the sum of the 
muon and antimuon events in each bin. Therefore, in this case we have just a 60 bin data. Such 
a non-magnetized iron detector could then rule out the fake octant solution at 3<r level for the 
normal mass ordering if sin 2 ^(true) < 0.372(> 0.625) for sin 2 6>i 3 (true) = 0.04 or sin 2 ^(true) < 
0.398(> 0.600) for sin 2 #i 3 (true) = 0.04. It is not surprising that the potential to rule out the 
fake octant deteriorates significantly in this case, since the matter effects crucial for the octant 
sensitivity get diluted if neutrino and antineutrino events were added together. Distinguishing 
the neutrinos from the antineutrinos is therefore extremely important for the physics potental of 
experiments sensitive to large matter effects. 

8 Discussions and conclusions 

The measurement of both the magnitude and sign of the deviation D of sin 2 #23 from its maximal 
value 0.5 is of utmost theoretical importance. The best current limit on this parameter D comes 
from the SK atmospheric neutrino experiment giving \D\ < 0.16 at the 3a level [5]. More precise 
measurements of #23 and hence of D are expected from future atmospheric neutrino data both from 
the currently running SK experiment and from the planned Megaton water Cerenkov detectors. 
Significantly better constraints are expected from data that will emerge from forthcoming long 
baseline experiments, owing to their larger statistics and lower systematic errors. However, in both 
these classes of experiments, 6*23 will be mainly determined by the (and/or v^) disappearance 
channel in vacuum, which predominantly depends on sin 2 2# 23 . This leads to two very important 
consequences. The first one relates to the fact that the mixing angle #23 is very close to being 
maximal. The fact that 5(sin 2 #23) equals <5(sin 2 2# 23 )/(4 cos 2^23) has the following implication. 
Even though one could determine the value of sin 2 2# 2 3 at the percentage level from the next 
generation long baseline experiments, the uncertainty in sin 2 #23 would still remain in the region 
of 10-20%, depending on sin 2 #23 (true). The second consequence of the predominant sin 2 26> 23 
dependence of the disappearance probability in long baseline experiments means that they are 
almost insensitive to the octant of #23 and hence to the sign of D. 

In this paper we have argued in a quantitative way that the presence of matter effects in the 
neutrino survival probability enhances its sensitivity to #23 . In particular, we have shown that 
the strong sin 2 #23 dependence of the matter effects in P w can be used to increase the sensitivity 
of z/^'s and P^'s disappearance experiments to D. Using the signal for atmospheric and in 
a large magnetized iron calorimetric detector like the planned ICAL at INO in India, we have 
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demonstrated the followoing fact: for large enough values of sin 2 #13 (true), the detection of matter 
effects in the data would be feasible and would lead to better constraints on \D\. Appropriate 
binning of the data in both energy and zenith angle holds the key to the observation of matter 
effects in the resultant signal. We have distributed our simulated events into twelve zenith angle 
bins covering both up and down going neutrinos - each of which was then further divide into five 
energy bins, giving us sixty bins of data. The charge discrimination capability of a magnetized 
detector like ICAL will allow the separation between v^- and z/^-induced events. Hence we have 
simulated separate data sets for z/^'s and z/ M 's. Since for a given neutrino mass ordering, substantial 
matter effects appear in only either the neutrino or the antineutrino channel, the effective use 
of separate data sets for and will enable the extraction of matter effects in the survival 
probability through a statistical analysis. We have defined a x 2 function to analyze the data set 
expected from atmospheric neutrino events collected in an ICAL-like detector and presented the 
C.L. limits on \D\. We have noted that the presence of matter effects increase the sensitivity of 
ICAL to \D\ somewhat. For sin 2 #13 (true) = 0.04, matter effects are seen to improve the 3a limit 
on \D\ from 0.090 to 0.084. We expect to measure \D\ within 18% (17%) at 3a with 1 MtonY 
data when sin 2 #13 (true) = 0.02 (0.04). This is comparable to the limit on \D\ expected from the 
forthcoming long baseline experiments and slightly better than what is expected from atmospheric 
neutrino experiments with Megaton water Cerenkov detectors. 

Matter effects in open up a new utility for an ICAL-like detector - sensitivity to the 
octant of 823- As is well known, the atmospheric neutrino data set collected so far by the SK 
experiment constrains the atmospheric neutrino mixing angle in the form sin 2 2#23 and is hence 
insensitive to the octant of #23 if its true value were non-maximal This ambiguity in whether #23 
is smaller or greater than ir/4 leads to an additional two- fold degeneracy [38] in the measurement 
of the mixing angle #13 and the CP phase 8 in long baseline experiments, looking for electron 
neutrino appearance in an original muon neutrino beam through the P^ e channel which depends 
on sin 2 #23 and hence to the octant of #23 • The P M/i channel in the long baseline experiments, which 
is expected to give the best limits on #23, is almost insensitive to its octant [15]. The proposed ways 
to tackle the #23 octant ambiguity in the determination of #13 and 5 in long baseline experiments 
include (1) combining the data from the P Me (or P efl ) channel of the long baseline experiment with 
data from the next generation reactor experiments [34], (2) combining the data from the P Me (or 
P e n) channel of the long baseline experiments with different energies and baseline lengths and (3) 
combined studies of the P Me or P eM (golden) and the P eT (silver) channels [39]. 

The octant ambiguity can be resolved from a direct measurement of the sign of D utilizing 
atmospheric neutrino data. It was shown in [9] that the atmospheric neutrino data in SK like 
experiments could be sensitive to the octant of #23 through the Am 2 x driven sub-dominant os- 
cillations. The latter are sizable for very small values of the neutrino energies and result in an 
excess in the sub-GeV electron sample in Megaton water Cerenkov detectors. They conclude 
that, with a statistics fifty times the current SK statistics, this effect can be used to determine 
the correct octant of #23 at the 3a level if the true value sin 2 ^23 (true) < 0.36 or > 0.62. We 
have shown in this paper that for large enough sin 2 #13 (true), the observation of matter effects in 
atmospheric neutrinos in an ICAL like detector could be used very efficiently to determine the 
octant of #23- The observation of significant matter effects in P w allows the rejection of the false 
octant solution of # 23 at the 3a level for sin 2 # 23 (true) < 0.40 (< 0.42) or sin 2 6> 23 (true) > 0.59 
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(> 0.57) for sin 2 0i 3 (true) = 0.02 (0.04) with 1 MtonY statistics in ICAL, given a normal neutrino 
mass ordering. For an inverted mass ordering, the false #23 can be rejected at the 3a level for 
sin 2 #23(true) < 0.37 (< 0.40) or sin 2 23 (true) > 0.62 (> 0.60) for sin 2 # 13 (true) = 0.02 (0.04). The 
sensitivity to sgn(D) in ICAL is therefore slightly less for an inverted mass ordering than for a 
normal one. But it is still better than what is expected from other experiments. Hence, using 
matter effects to pick the right octant in an INO-like experiment seems to be the most promising 
way of resolving the #23 octant ambiguity by determining the sign of D. 

We thank S. R. Dugad, S. Goswami, M. V. N. Murthy, D. P. Roy and A. Raychaudhuri for their 
helpful comments. 
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Figure 1: as a function of L with fixed E (and at various values of E) for neutrinos travelling 
in vacuum and matter. The black solid lines and the dashed green lines are for propagation in 
matter with sin 2 # 2 3 = 0.5 and 0.36 respectively. The red dotted lines and the blue dot-dashed 
lines are for that in vacuum for the same respective values of sin 2 #23 . The vertical dashed lines 
represent the mantle-core boundary inside the earth. 
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Figure 2: The average E res (blue solid line) as a function of the baseline length in earth. Also 
shown are ^spmini (red dotted line), -Espmax (green dashed line) and -Espmin2 (black dot-dashed 
line) in vacuum. They curve owing to the logarithmic scale of the horizontal axis. 
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Figure 3: Plots of -P^ 2,3 , the notional quantities of Eq. (27) defined in the limit a — > 0, as functions 
of the neutrino energy E with a fixed baseline length L for different L- values as shown. The black 
solid lines, red long-dashed lines and green dashed lines correspond to u^s travelling in matter 
with sin 2 # 2 3 = 0.5, 0.4 and 0.6 respectively and the blue dotted lines to the same in vacuum with 
sin 2 #23 = 0.5. 
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Figure 4: as a function of E for fixed values of L. The black solid and the dashed green lines 
are for neutrinos travelling in matter with sin 2 #23 = 0.5 and sin 2 #23 = 0.36 respectively. The red 
dotted and the blue dot-dashed lines show the corresponding respective cases in vacuum. 
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Figure 5: The up-down asymmetry expected for muon neutrinos in energy bins of width 2 GeV. 
The six panels show the data for six different zenith angle (or L) bins corresponding to upward 
neutrinos travelling between L m = - 2000 km, L m = 2000 - 4000 km, L m = 4000 - 6000 km, 
L m = 6000 - 8000 km, L m = 8000 - 10000 km and L m = 10000 - 12000 km respectively inside the 
earth, cf. Table 3. The solid black lines and the solid magenta lines are for neutrinos travelling in 
matter with sin 2 623 = 0.5 and 0.36 respectively. The dashed black lines and the dashed magenta 
lines are for neutrinos travelling in vacuum with sin 2 6*23 = 0.5 and 0.36 respectively. For all cases 
we have taken Am| x = 2 x 10 -3 eV 2 and the benchmark values of Table 2 for the other parameters. 
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Figure 6: The difference between the up-down ratio for the neutrinos (Un/Dn) and antineutri- 
nos (Ua/Da) shown for the various energy and zenith angle bins. The solid black lines are for 
neutrinos/antineutrinos travelling in matter with sin 2 # 2 3 = 0.5, while the solid magenta lines are 
neutrinos/antineutrinos travelling in matter with sin 2 623 = 0.36. The dot-dashed blue line shows 
Un/Dn — Ua/Da = for reference. The other oscillation parameters are chosen as in Fig. 5. 
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Figure 7: The regions of Amg 1 (true) and sin 2 # 23 (true) where maximal 6 2 z mixing can be rejected 
by using 1 MtonY (upper panels) and 3.37 MtonY (lower panels) atmospheric neutrino data in 
ICAL at la (white band), 2a (blue band) and 3a (green band). The hollow dark lines show the 
corresponding bands for neutrinos travelling in pure vacuum. Benchmark parametric values of 
Table 2 have been assumed. 



34 




-1 -0.5 

cos z 



Figure 8: The zenith angle distribution of muon events for four energy bins shown in the four 
panels and for Am^ = 2 x 10~ 3 eV 2 . The black solid lines show the events for maximal sin 2 #23, 
the dotted magenta lines show the events for sin 2 #23 = 0.4 and the blue dot-dashed lines show the 
events for the unoscillated flux. 
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Figure 9: Same as Fig. 8, but for Am 2 x = 4 x 10 3 eV 2 . 
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Figure 10: \ 2 as a function of sin 2 # 2 3, showing the "octant sensitivity" of ICAL. The three panels 
project the expected sensitivity of ICAL for three different true value of sin 2 6*23: sin 2 6*23 (true) = 
0.42 (left panel), sin 2 #23 (true) = 0.46 (middle panel) and sin 2 2 3 (true) = 0.54 (right panel). The 
magenta long-dashed lines, blue short-dashed lines, green dot-dashed lines, red dotted lines and 
black solid lines are for sin 2 # 13 (true) = 0.00, 0.01, 0.02, 0.03 and 0.04 respectively. 
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Figure 11: A% 2 as a function of sin 2 #23 (true), showing the octant sensitivity of ICAL for the 
normal (left panel) and inverted (right panel) neutrino mass ordering. The blue short-dashed 
lines, green dot-dashed lines, red dotted lines and black solid lines are for sin 2 6*13 (true) = 0.01, 
0.02, 0.03 and 0.04 respectively. 
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